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Abstract — A construction of a new family of distributed space 
time codes (DSTCs) having full diversity and low Maximum 
Likelihood (ML) decoding complexity is provided for the two 
phase based cooperative diversity protocols of Jing-Hassibi and 
the recently proposed Generalized Non-orthogonal Amplify and 
Forward (GNAF) protocol of Rajan et al. The salient feature of 
the proposed DSTCs is that they satisfy the extra constraints 
imposed by the protocols and are also four-group ML decodable 
which leads to significant reduction in ML decoding complexity 
compared to all existing DSTC constructions. Moreover these 
codes have uniform distribution of power among the relays as 
well as in time. Also, simulations results indicate that these codes 
perform better in comparison with the only known DSTC with 
the same rate and decoding complexity, namely the Coordinate 
Interleaved Orthogonal Design (CIOD). Furthermore, they per- 
form very close to DSTCs from field extensions which have same 
rate but higher decoding complexity. 

I. Introduction 

Cooperative diversity is a technique by which multiple 
terminals (users or relays) cooperate to form a virtual antenna 
array thereby leveraging the spatial diversity benefits even if a 
local antenna array is not available. A cooperative diversity 
protocol dictates how the users would actually cooperate 
among themselves to achieve the required diversity order. 
Several cooperative diversity protocols have been proposed in 
the literature [l]-[6]. In this paper, we focus on the two phase 
based protocols of Jing-Hassibi [1] and the GNAF protocol 
of Rajan et al [2] for three reasons- (i) the operations at 
the relay nodes are considerably simplified, (ii) we can avoid 
imposing bottlenecks on the rate by not requiring the relay 
nodes to decode and (iii) the framework of distributed space- 
time codes allows for more flexibility and higher spectral 
efficiency [1], [2], [4]. Transmission in this protocol comprises 
of two phases- broadcast phase and cooperation phase. In the 
broadcast phase, the source broadcasts its information to the 
relays and the destination. In the cooperation phase, each relay 
transmits a linearly processed version of the received vector. 
For this purpose, each relay is equipped with a unitary matrix 
which we call 'relay matrix'. It was shown in [1], [2] that 
to the destination it would appear as if a space-time code 
was transmitted from colocated multiple antennas. Further the 
design criteria to achieve full diversity also remains to be the 
well known rank criteria for colocated MIMO (Multiple Input 
Mutiple Output Systems). However it is important to note that 
DSTCs need to satisfy many additional constraints on the 
code structure (for example, unitary relay matrices) and are 



different from the traditional Space-time codes for colocated 
MIMO systems, which are designed without respecting such 
constraints. 

The ML decoding of a STBC in K complex variables 
Xi,X2,--- ,xk is, in general, joint decoding of all the K 
variables. However, for the Alamouti code K = 2 and the 
variables x\ and X2 can be decoded independently for ML 
decoding. In general, if K = gX and the variables can be par- 
titioned into g subsets each containing A number of variables 
and the ML decoding can be done for the variables of a subset 
independently of the variables of other subsets the code is said 
to be g-group ML decodable or A-symbol decodable [9], [10]. 
The Alamouti code is single-symbol decodable or two-group 
ML decodable. Following the work of [1], several distributed 
space time codes [6], [7], [8] were proposed. However, most of 
these code constructions did not consider the important aspect 
of ML decoding complexity at the destination. This problem 
gains significant importance especially if the number of relays 
in the network is large. An initiative in this direction was 
first taken in [11] wherein two-group ML decodable DSTCs 
were proposed using a construction procedure called 'doubling 
construction'. Later in [3], a class of rate one, full diversity, 
four-group ML decodable DSTCs called Precoded CIODs(Co- 
ordinate Interleaved Orthogonal Designs) were proposed for 
arbitrary number of relays. However these DSTCs had a large 
number of zero entries in the design which led to a large Peak 
to Average Power Ratio (PAPR). In [10], a class of four-group 
ML decodable STCs were proposed for the colocated MIMO 
systems. However these STCs fail to satisfy the additional 
constaints on the code structure imposed by the cooperative 
diversity protocol. Hence it is important to address the problem 
of constructing a new class of four-group ML decodable 
DSTCs which have low PAPR and uniform distribution of 
power among the relays. 

The main contribution of this paper is in providing a new 
construction of DSTCs with the following salient features. 

• Source can transmit \ complex symbols per channel use 

• Full diversity 

• Four group ML decodable 

> Unitary relay matrices. (This eliminates the need for a 
whitening filter at the destination thus further reducing 
decoding complexity.) 

« Unitary weight matrices together with unitary relay ma- 
trices makes the power distribution uniform among the 



relays and in time and hence results in low PAPR. 

The paper is organized as follows: We briefly describe the 
notion of four-group ML decodable codes and also state the 
DSTC design constraints which are due to the protocol in 
Section [TT] In Section [ill] we describe the code construction 
procedure and show that they possess the salient features 
stated above. Few illustrative examples of DSTC construction 
are also provided for 4 and 16 relays. Simulation results 
are presented in Section [IV] Conclusions and discussions on 
further work comprise Section [V] 

Notation: For a complex matrix A, A*, A T and A H denote 
the conjugate, transpose and conjugate transpose respectively. 
Aj denotes the real matrix obtained by taking the real parts of 
all the entries of the matrix A and Aq denotes the real matrix 
obtained by taking the imaginary parts of all the entries of the 
matrix A. For a square matrix B, \B\ and Tr (B) denote the 
determinant and trace of the matrix B respectively. 

II. Four group ML decodable DSTC design 

PROBLEM 

In this section, we briefly describe the two phase based 
cooperative diversity protocol of [1], [2] and introduce the 
problem statement. Consider a wireless relay network consist- 
ing of a source node, a destination node and R other relay 
nodes which aid the source in communicating information to 
the destination. The channel path gains from the source to the 
i th relay, denoted by fi and those from the j th relay to the 
destination denoted by gj are all assumed to be i.i.d CAf(0, 1). 
The channel path gain, go from the source to the destination is 
also assumed to be CAf(0, 1). All the nodes are equipped with 
a single antenna and are subject to the half-duplex constraint, 
i.e., a node cannot transmit and receive simultaneously. Fur- 
ther, we assume that the nodes are synchronized at the symbol 
level. Each transmission from source to destination comprises 
of two phases- broadcast phase and cooperation phase. In 
the broadcast phase, the source transmits a T length vector 
s taken from a codebook consisting of information vectors 
c i = {si, . . . , sl} satisfying E {s ff s} = 1 to all the relays 
and the destination. In the cooperation phase, all the relay 
nodes are scheduled to transmit together (assuming symbol 
level synchronization) a DSTC. For this purpose, each relay 
is equipped with a unitary matrix A; which we call 'relay 
matrix' . To be precise, the i th relay transmits a scaled version 
(to satisfy power constraint) of A^i or A%r\, where denotes 
the received vector at the relay. The signal model is as shown 
in ([T]i at the top of the next page, where 

> ti denotes the vector transmitted by the i th relay and i>i 
denotes the additive noise at the relay whose entries are 
assumed to be i.i.d CAf(0, 1). 

- ud.i and yn,2 denote the received vector at the desti- 
nation during the broadcast phase and cooperation phase 
respectively. w\ and W2 represent the additive noise at 
the destination whose entries are i.i.d CAf(0, 1). The 
quantities ir\ and TT2 are the power allocation factors 
satisfying tti + ^Ti? = 2T so that P represents the 



total average power spent by the source and the relays 
together. 

The received vector at the destination can be written in 
matrix form as follows 



Vd,i 
Vd,2 



7T27T1P 2 



SH + W 



(2) 



V 

where S, H and W are as shown in (O at the top of the 
next page. The DSTC in this case is the collection of all the 
(2T) x (P+l) matrices S. Observing the structure of S in (0, 
we see that it is sufficient to design the submatrix of S given 
by S E r =[ A lS ... A N s A N+1 s* . . . A R s* ] . We 
are interested in the case when the submatrix Ser is obtained 
from a linear dispersion STBC say S(X). Let S(X) = 
Z^iLi XiAi where, x\, X2, ■ ■ ■ , xk are the K real variables of 
the linear STBC S(X) and the matrices A, e C TxNt , called 
the 'weight matrices', define the code. 

The vector X = [xi,x%, . . . , xk] T G srf C 
information symbol vector. 

Suppose we partition the 
tor as X T = [Xf,X. 

fo fe +i>*3 fe +2, • • • ,%j k +n k ], ji = and j k = J2i=i n i for 
k = 1, 2, . . . , g and their corresponding set of weight matrices 
also into g groups L k , k = 1, ... ,g, the k th group containing 
n k matrices, then S(X) can be written as 

a 

S{X) = Y J S k {X k )^CTC S k {X k ) 



i K is called the 



information symbol vec- 
2>--> X g]> where 



k=l 



1=1 



i.Aj 



Now if the information symbols in each group take values 
independent of information symbols in the other groups and 
if the weight matrices satisfy 



AfAj + AfAi = 0, Vi E L p ,Vj eL„p^ q. 



(4) 



then, it can be easily shown [9], [10] that S(X) is (/-group ML 
decodable. In other words, the ML decoding can be performed 
by minimizing the metric 



II V 
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P 2 



TTlP+1 



Sk(X k )H 



(5) 



for each 1 < k < g individually instead of minimizing the 
computationally more intensive metric 



D 



'TTaTTlP 2 
7TlP+ 1 



S{X)H 



(6) 



Hence the ML decoding complexity is reduced to a large 
extent depending on the value of g. In particular, the decoding 
complexity is reduced from 2 A to g2~a where, A depends 
on the rate of transmission as measured in bits/s/Hz. In this 
paper, we consider only the case of g = 4. Combining all the 
given requirements, the 4-group ML decodable DSTC design 
problem is thus to find space-time block codes satisfying the 
following three constraints. 

1) Any column should contain only the variables of the 
design or only their conjugates. 



Vd,x = V^Pg Q s + wi 

r% = y/ niPfjS + Vj,V i = l,...,R 

t< = v^SS 4 ^' i=1 '---> N (i) 

VD,2 = Yli=l 9%U + w 2 
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9n!n gN+if^+i 

Wi 



9rJe 



(3) 



ttTp+I (Sill 9iAiVi + X^ilw+l 



w 2 



2) All the relay matrices should be unitary. 

3) The weight matrices of the code should satisfy the 
conditions for 4-group decodability. 

III. Code Construction procedure 

Having described the problem statement, in this section we 
explicitly construct a new class of rate one, full diversity four- 
group ML decodable DSTCs. 
Consider the following design of size R x R, 

' A -B H ' 
B A H 



(7) 



where, A and B are identical | X | designs in different 
variables. If the codewords in A commute with the codewords 



in B, then we have 



S H S 



A H A + B H B 

BB H + AA H 



(8) 



This ensures that the code S is two-group ML decodable, one 
group involving variables in A and the other group involving 
the variables in B. This fact was exploited in [11] to construct 
two-group ML decodable DSTCs. This construction procedure 
was named as 'Doubling construction' in [11]. However, note 
that in addition if A and B are individually two-group ML 
decodable, then it can be easily shown that the code S will be 
four-group ML decodable. In this paper, we essentially exploit 
this fact. Hence we require two-group ML decodable codes 
whose codewords commute among themselves. Towards that 
end, consider the following 2x2 designs. 



si s 2 
-s 2 Sl 



and C2 



si s 2 
s 2 si 



We can show that both C\ and C 2 are two-group ML decod- 
able codes. For the design C 2 , the two groups are {su, s 2 q} 
and {siq,s 2 j} respectively. More importantly, we have the 
following identity true for any complex number 7. 



Sl 


s 2 




Sj s' 2 




Sj s' 2 




Sl 


S2 


_ 7S2 


Sl 




. 7S2 s[ _ 




_ 7S 2 si _ 




_ 7S2 


Sl 



In other words, codewords arising from design C± or C 2 
commute. Thus we can construct a 4 x 4 design by substituting 
either C\ or C 2 for A and B in (0 as shown below. 



S 



si 

S2 
S3 
S4 



S2 
Sl 
S4 
S3 



(9) 



It can be easily verified that S is indeed 4-group ML decod- 
able. The four-groups are {su,s 2I }, {si Q ,s 2Q }, {s 3I ,s 4I } 
and {S3Q, S4q} respectively. To extend this approach for more 
number of relays, we need to find designs analogous to C\ 
and C 2 for higher dimensions. Towards that end, consider the 



Y x design 



D 



W 
X 



X 
W 



where, W and X are identical | x | designs in different 
variables. We call this construction as the 'ABBA' construc- 
tion. It is easy to check that if the design W or X is g- 
group ML decodable, then the design D obtained using ABBA 
construction is also g-group ML decodable. Further, we have 
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X 
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WX' + XW 1 
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Thus we observe that codewords obtained from design D 
commute if WW = WW, X'X = XX', WX = XW 
and X'W — WX'. Since W and X are taken to be identical 
designs, this simply means that if codewords in W commute, 
then codewords of D also commute. Using the above facts, we 
now give our construction procedure for any R = 2 a relays 
for a > 2. 

Construction 1: The steps in the construction procedure are 
enlisted as below. 

Step 1: Starting with either C\ or G 2 , keep applying ABBA 
construction iteratively on it till a|x| design C is obtained. 
The design C will be two-group ML decodable and it will also 
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have commuting codewords. 

Step 2: Apply doubling construction (O to C by substituting 
C for A and £>(albeit with different variables). Thus we obtain 
a Rx R design which will be four-group ML decodable. 

Note that the codes designed using Construction Q] have the 
property that any column has only the variables or only their 
conjugates. Moreover it can be shown that these codes have 
unitary relay matrices as well as unitary weight matrices. This 
leads to uniform distribution of power among the relays as well 
as in time and hence low PAPR. 

Example 1: Using construction [T] we get the 
16 x 16 design shown in ( TTOb at the top of this page. 
The four groups of variables that can be decoded 
separately are {s u , s 2 j, . . . , s 8 j}, {s 1Q , s 2 q, ■ ■ ■ , s 8 q}, 
{sgi, sio/, ■ ■ ■ , siei} and {s 9Q , s 10 q, . . . , s 16Q } respectively. 
Observe that any column has only the variables or only its 
conjugates. Further note that all the relay matrices are unitary. 
It can also be checked that all the weight matrices are unitary. 

A. Diversity of constructed codes 

The constructed codes can be made to achieve full diversity 
by letting the variables in each of the four groups to take 
values from an appropriately rotated Z~ lattice constellation. 
We would like to emphasize here that the variables in a group 
should be allowed to take values independent of the variables 
in the other groups. Algebraic number theory provides effec- 
tive means to construct rotated Z™ lattices with full diversity 
and large minimum product distance [12]. Such lattices can 
be used to obtain full diversity for our codes. We will now 
illustrate how this can be done for the 4x4 design shown in 
©. 

Note that the matrix S is invertible if either A or B is 
invertible. This is because, 

|5| = \S H S^ = \{A H A + B H B){AA H + BB H )\i 
> m&x{\A\ 2 ,\B\ 2 } 

(ID 

where, the second inequality was obtained on application of 
Corollary 4.3.3 in [13]. Therefore it is sufficient to show how 
Ci can be made fully diverse. We have 



Therefore the condition for full-diversity is Asi ^ ±As2. 
Thus, if Asi/ 7^ ±As2/ and Asiq ^ ±As2Q then full 
diversity is guaranteed. Let us define p = -A=(s ir + s 2 /) and 
q = -jf=(su — S2i). Then the required condition is simply 
ApAq 0, which is nothing but non zero product distance. 
So let p, q take values from a rotated Z 2 lattice constellation 
designed to maximize product distance [12]. Then let su = 
and s-2i = 2=£. Similarly it is done for siq,S2q also. 
Thus we have 
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where, G 



\AC 2 \ = (Asi + As 2 )(Asi - As 2 ). 



(12) 



-0.5257311121 -0.8506508083 
_ -0.8506508083 0.5257311121 
In this manner, we can show that it is possible to achieve 
full diversity for all the codes in this paper. 

IV. Simulation results 

In this section, we compare the performance of the newly 
proposed DSTC for a 4 relay network (0 with that of DSTCs 
from CIOD [2] and field extensions [6], [7], [8]. The power 
allocation factors used for the simulation are 7Ti = T, and 
7r 2 = Jj. The signal sets chosen for the different DSTCs are 
as follows: 

1) CIOD - QPSK rotated by 31.7175° 

2) Field extension - QPSK 

3) Newly proposed code - QPSK appropriately rotated as 
explained in subsection A of Section Hill 

FigQ] shows the codeword error rate comparison of the 
newly proposed DSTC with that of DSTCs from CIOD and 
field extensions. Observe from FigQ] that the newly proposed 
code performs slightly better compared to the 4-group ML 
decodable CIOD [2]. Moreover its performance is only slightly 
inferior compared to the 1 -group ML decodable DSTC from 
field extensions [6], [7], [8]. 

V. Discussion 

A class of full diversity 4-group ML decodable DSTCs 
with unitary relay matrices and unitary weight matrices were 
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Fig. 1 . Performance comparison of the proposed code with DSTCs from CIOD and field extension 



explicity constructed in this paper for number of relays of 
the form R = 2 a ,a € N. The effect of power allocation on 
the error performance was not considered in this work and 
is an interesting problem for further work. Further study on 
the codes of this paper has revealed an interesting algebraic 
structure which gives a unifying theory of several known codes 
in the literature including the ones in this paper. An algebraic 
approach to designing these codes is currently under progress. 
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